Effect of fractal disorder on static friction in the Tomlinson model 
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We propose a modified version of the Tomlinson model for static friction between two chains of 
beads. We introduce disorder in terms of vacancies in the chain, and distribute the remaining beads 
in a scale invariant way. For this we utilize a generalized random Cantor set. We relate the static 
friction force, to the overlap distribution of the chains, and discuss how the distribution of the static 
friction force depends on the distribution of the remaining beads. For the random Cantor set we 
find a scaled distribution which is independent on the generation of the set. 
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I. INTRODUCTION 

Friction is a long studied phenomenon in physics [H-Q- 
Still, the regime of validity and the microscopic origin 
of the empirical laws of static and dynamic friction (that 
goes under the name of Amontons-Coulomb) are not well 
established. But the advancements of technology in the 
last few decades has triggered both theoretical jj-B, as 
well as experimental investigations in this field |10L . 

Surfaces which appear smooth may contain roughness 
at the micrometer scale [r2| as well as impurities in the 
contact region [l3| which will influence the friction prop- 
erties of the materials. If we limit ourselves to the study 
of atomically clean surfaces, it is sufficient to consider 
the relative motion of the atomic layers in contact. The 
earliest attempt to model such a situation was carried 
out by Tomlinson [Til- Tomlinson's model is a chain 
of beads, representing atoms, all of which are individu- 
ally attached to a body above. The chain is dragged on 
a periodic potential representing a corrugated substrate 
(see Fig. Q] (a)). Each bead interacts with the potential, 
but there is no interaction among the beads. As in the 
Tomlinson model for dry friction (see e.g. [Hj]), the sur- 
face atoms are considered to be independent oscillators 
capable of absorbing a finite energy and momentum. Al- 
though they are connected by the chain, the energy is 
assumed to get dissipated within the bulk and not trans- 
ferred to the neighbouring atoms along the chain. This 
model has regained its importance in recent years be- 
cause, in experiments using atomic force microscopy, the 
contact probe can be studied at an atomic level (see e.g., 
0]). There are also a number of similar models using 
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chains of beads moving over a substrate. Most notably is 
the Frenkel-Kontorova model [l6jj, where the beads also 
interact with the nearest neighbors. 

In this paper, we propose a modified version of the 
Tomlinson model. We consider a two chain version of 
the model, where one of the chains slides on top of the 
other (see Fig. [l](b)). We study the effect of disorder 
in the static friction force. In particular, we introduce 
vacancy disorder, by removing beads, in both chains. We 
relate the static friction force to the measure of overlap 
of beads. The disorder in the chains is such that the 
remaining beads are distributed in a scale invariant way. 
For this we use a generalized version of the Cantor set. 
This leads to a certain kind of overlap distribution, and 
we compare this to other ways of distributing the beads, 
and discuss how this relates to the distribution of static 
friction force. 

It may be mentioned at this point that rough surfaces 
have a self-affine character and therefore can be repre- 
sented by self-similar fractals [13] ■ We represent the 
disorder in the Tomlinson model by a generalized ran- 
dom Cantor set, which is one of the simplest examples 
of a fractal. Although the Cantor set is a purely mathe- 
matical construction, it has been studied extensively by 
physicist as it encapsulates the essential property of self- 
similarity and scale invariance. 

The paper is organized in the following way: In section 
(II) we define the model. Then in section (III) we briefly 
discuss the basic definitions and useful properties of the 
Cantor set, which will be required in the subsequent dis- 
cussions. In section (IV) we find the expressions for the 
static friction force distribution for Tomlinson's model 
with beads distributed as random cantor sets. We con- 
clude by summarizing the results in section (V). 
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II. THE MODEL 

The static friction force F s , is defined to be the mini- 
mum force needed to initiate sliding between two objects 
in contact. If the applied force exceeds this threshold 
F > F s , the objects will move relative to each other. 
The presence of this force is due to a local minimum 
in the energy landscape of spatial configurations of the 
two objects. Two disordered surfaces may however have 
several configurations which correspond to local minima, 
and the static friction may have different values in the 
different configurations. 

We consider a two chain version of Tomlinson's model, 
where one chain slides on top of the other (see Fig. [T] 
(b)). We introduce defects in terms of vacancies (re- 
moved beads). The vacancies are introduced randomly, 
and are statistically identical in both chains. The spac- 
ings between the sites in the chain (either occupied by a 
bead or removed), are constant and equal in both chains. 

The underlying chain gives rise to a substrate potential 
for the chain above. We assume that the interaction po- 
tential between the beads in the opposite chains are short 
ranged and attractive, such that the substrate potential 
is on the form of a series of potential wells, correspond- 
ing to the remaining beads of the chain below (see Fig. 
□ (c)). 

The potential energy of a configuration of this modi- 
fied Tomlinson model will be determined by the sum of 
beads locked in potential wells. The static friction force 
of one bead locked in a potential well is determined by 
the maximum value of the derivative of the potential. We 
will not specify the functional form of the potential well, 
but we assume that the static friction force of one bead 
in a potential has a given value. The static friction force 
of a chain will therefore be directly proportional to the 
number of beads locked in potential wells. 

If there were no disorder, the overlap would be con- 
stant and the static friction force per bead would also 
be constant. If, however, there is disorder in terms of 
vacancies, the number of locked beads (referred to as the 
overlap) will not be constant but will have a set of pos- 
sible values. These values, and hence the static friction 
force, will follow a distribution. 

We assume that the expected number of overlapping 
beads is given and investigate how the distribution vary 
around the expectation value. We will see that this dis- 
tribution will have certain features if the beads (and the 
potential wells from the substrate) are distributed with 
a spatial correlation between them. This distribution 
will differ significantly from the case when the remaining 
beads (and the potential wells from the substrate) are 
distributed uniformly. 

We are especially interested in the properties of the dis- 
tribution when the beads arc distributed with a spatial 
correlation which decays as a power-law in the distance 
between the beads, as this correspond to a fractal disor- 
der in the chain. A binary chain (with vacant or occupied 
sites) is, however, not uniquely determined by the spatial 



correlation structure. We choose to utilize generalized 
random Cantor sets for the distribution of the remaining 
beads (see Fig. [3]), as the random Cantor set has the 
desired correlation structure on an average (this is a re- 
sult of the self-similar property of the set). Moreover the 
structure of the random Cantor set is simple enough to 
allow theoretical results for the overlap distribution. 
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FIG. 1: (a) The Tomlinson model: A chain of beads 
being dragged with an applied force F in a periodic 
potential V(x). (b) a two chain version of the 
Tomlinson model, with defects in terms of vacancies, 
(c) Same as (b) with the potential from the underlying 
chain illustrated. 



III. CANTOR SETS 

The prototype example of a fractal is the Cantor set C. 
The Cantor set is constructed by first removing the mid- 
dle third of the base interval [0, 1]. From each of the 
remaining intervals, [0, 1/3] and [2/3, 1], a middle third 
is again removed. The process of removing middle thirds 
of the remaining intervals is continued ad infinitum. The 
intervals which are left after the middle third of every re- 
maining interval have been removed n times, is referred 
to as the Cantor set C n , at generation n. C n becomes a 
true fractal as n goes to infinity C = lim C„. 

There are two self-similar transformations related to 
the Cantor set C. The transformations are (x) — x/3 
and SW(x) = (x + 2)/3. We can describe the process of 
removing the middle third of every interval by the action 
of these transformation, 
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C n = S^\C n - l )\JS^>(C n - 1 ) 

And wc can define C to be the subset of [0, 1] which is FI G. 3: Construction of a random Cantor set, s = 5 and 
invariant under the union of these two transformations. r = 3. 
Fig. [2] illustrates the construction procedure. 





FIG. 2: Construction of a triadic Cantor set. 



In order to investigate properties which emerge from 
the fractal nature of a subset of an interval, we will extend 
the notion of a Cantor set. We will construct a general- 
ized Cantor set by the action of r self-similar transforma- 
tions on the interval [0, 1]. Let a = {ai,<22, ...,a r } be a 
set of integers such that 1 < r < s and ttj < a 3 for i < j. 
The self-similar transformations be on the following form: 



And we will obtain the generalized Cantor set in the ob- 
vious way 

C„ = UsW(C„_x). 

i 

The Cantor set constructed this way will consist of r n 
line elements of equal length that may or may not be 
connected. We identify the integers in a as the positions 
of the remaining line elements at the first generation. 
The Box Counting dimension of this generalized Cantor 
set is given by Db = log(r)/ log(s) (see e.g. [13): an d 
the regular triadic Cantor set is obtained by taking s = 3 
and a = {1, 3}. 

The notion of a random Cantor set is ambiguous as 
we can randomize in several ways. We will define the 
random version of this generalized Cantor set analo gou s 
to Falconer's definition of the triadic Cantor set in [la |. 
For each generational step, divide every remaining line 
element from the previous generation into s elements of 
equal length and remove all but r of them. For each line 
segment, and for each generation, randomize the posi- 
tion of the r remaining intervals, but keep the number of 
remaining intervals fixed. The self-similar properties of 
this set hold only in a statistical sense. 

The generalized Cantor set will at generation n con- 
sist of r n of the s" equal distant parts of the inter- 
val [0,1]. This enables us to represent it as a vector 
(C n ,i,C n ,2, .••|C B , J »). The vector element C n ^ represent 
the ith interval of the Cantor set such that C n ,i takes 
the value 1 if the interval [(i — l)/s,i/s] is contained in 
C n , and the value if it is removed. The vector contains 



therefore s n elements, out of which r n takes the value 1, 
and the rest takes the value 0. 

The Cantor set is of measure zero. That is, the size 
of the Cantor set (sum of the remaining line segments), 
will go to zero as the generation goes to infinity. For that 
reason wc study the Cantor set at a finite generation, at 
the atomic level there is, of course, no way to continue 
the removal procedure. 



IV. OVERLAP 

A. Overlap of random sets 

Let us first consider the overlap of two chains with ran- 
domly placed vacancies, uniformly distributed along the 
chain. The overlap distribution is trivial, but we include 
it here as it will be instructive to compare it with the 
other results. 

Consider two independently generated vectors 
(Xi, ...,Xjv), (X[, ...,X' N ) containing N binary variates 
representing the chain (0 represents a removed bead 
at a given site). Each element takes the value 1 with 
probability p. The probability of an overlap at element i 
is therefore Pr(A"i = 1, X[ = 1) = p 2 . The overlap of the 
two vectors is given by Y = J^XiXj. And the overlap 
distribution is on the form of a binomial distribution: 



Pr{Y = x)=[ )p 2x (l-p 2 ) N - x . (1) 



B. Overlap of random Cantor sets 

In order to describe the overlap distribution for two ran- 
dom (independently generated) Cantor sets, C n ^ and 
C' n j, it is instructive to consider first the overlap of two 
random triadic Cantor sets (s = 3, r = 2), without a rel- 
ative displacement as illustrated in Fig. H](a). Notice, 
however, that Fig. [¥] also illustrates the construction of 
the Cantor sets. The figure does not suggest that the 
Cantor set is embedded in two dimensions. The over- 
lap at generation n is given by Y n = ~Y^C n ,iC' n i . For 
n = 1 one can easily calculate, by summing up the pos- 
sible outcomes, that the probability distribution of the 
overlap takes the following form 
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FIG. 4: (Color online) The overlap of two 
(independently generated) random Cantor sets at 
generation 3 (s = 3, r = 2). The generation number of 
the two sets are marked on the right side, (a) shows the 
overlap when the sets are perfectly aligned, and (b) 
shows the overlap when the sets are shifted relatively by 
an offset q. The overlap size is the number of 
overlapping elements at the final generation of the 
opposite sets (in the shaded region). The number of 
overlapping elements are 3 for both (a) and (b) . 



Call A the event that only one element overlap at the 
first generation of the construction, and B the event that 
two elements overlap at the first generation of the con- 
struction. A and B are mutually exclusive events, and 
P(A) + P(B) = 1. Consider the case when {C n ,i) is con- 
structed by an iterative procedure like in Fig. [3J We can 
write the following relation 

Pr(F„ =x) = \ ?<Yn =x\A) + ^ Pr(Y n = x\B). 
Note that 

Pr(F„ = x\A) = Pr(Y„_ x = x), 

Pr(F„ = x\B) = ]T Pr(y rl _! =x-i) Pr(Y„_ x = i). 

So we have 

Pr(F„ =x) = 

I Pr(y„_! = x) + 1 ]T Pr(r„_! = x - i) Pr(y„_ x = i). 



As n grows large, we can approximate this by a con- 
tinuous distribution f n {x) = Pr(K„ = x). The equation 
above translate to the following equation: 

/„(x) = \f n -l{x) + * fn-l)(x), (2) 

where (a * b)(x) is the convolution a(y)b(x — y) dy, 
and f n (x) = for x < 0. The overlap difference be- 
tween two consecutive generations in the limit of large 
n, should not alter the qualitative behavior of the prob- 
ability distribution. We therefore assume that there is a 
n-independent distribution / such that 

f n (x) = f(x/fi n )/Li n , (3) 

for a scaling variable (i. This is verified by numerical 
calculations, with /j, = (Yi) = r 2 /s = 4/3 which is the 
expected overlap at the first generation. f(x) should then 
obey the equation 

/(3/AO/M =§/(*) + £(/ */)(*)■ ( 4 ) 

To proceed with this functional equation we can take 
the Fourier transform on both sides. Using the elemen- 
tary properties of the Fourier transformation we get 

If we assume that f(k) is analytical in the complex 
plane, we can expand it as f(k) =Y1 ■ By matching 
terms we get a family of solutions 

^ n— 1 

a = 1, ax € C, a n = 4 ^ 7I _ 1 _ ^ % a »-i> n > 2 - 

The real part of f(k) has to be even, and the imaginary 
part has to be odd to make f(x) real. We therefore choose 
ai to be purely imaginary. 

Note that the Fourier transform is identical to the 
characteristic function tp(k) = (exp(ifca;)) of f(x). The 
moments of the distribution is given by the relation 
(x m ) = (-i)"V m )(0). The first moment of f(x) is by 
construction one, (x) = 1, so we choose a\ = i. (do = 1 
is consistent with the normalization of f(x).) All other 
moments are specified by 

(x m ) = (-i) m m\a m , (5) 

and the probability distribution /(x), is uniquely deter- 
mined in terms of its moments. 

For the general random Cantor set, defined in the pre- 
vious section, we have the following results. For the over- 
lap probability for the first generation we have 
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Pr(y 1 (s ' r) = x) =p (s ' r) (x) 



cJ \i — xJ 



Analogous to the distribution in eqn. ([3]), we have that 
the overlap distribution for the general case should be of 
the form 



3=1 



(*), (6) 



V 



j — 1 times 

and the n invariant distribution f( 8,r J(x) is given by 



/ 



J2p 

3 = 1 



(s,r) 



(./) 



p(s,r) 



*/ 



( s ,r) 



\ 



j — 1 times 



(7) 

with \i = r 2 1 s. By expanding the Fourier transform of 
j( s ' r )(x) as an analytical function, we can specify the 
moments for the probability distribution as we did in 
eqn. ([5]), but the values of s and r have to be specified 
to do so. 

Up until now we have considered the overlap of the ran- 
dom Cantor sets without any offset, i.e. the Cantor sets 
are placed such that the first element of the top Cantor 
set (vector element C ni i) overlap the first element of the 
lower Cantor set (see Fig. U (a)). In order to study the 
overlap when we also have a random offset, we assume 
periodic boundary conditions as indicated by Fig. |U(b). 
There is no obvious way to construct an analytical ex- 
pression for the distribution with the randomized offset, 
but we can generate the distributions numerically. Fig. 
[5] shows how the overlap distribution converges to the 
n-invariant distribution with and without a randomized 
offset, for (a) s, r = 3, 2, and (b) s, r = 6, 5. 

Fig. [5] (a) shows the distributions for different values 
of s with r = s — 1, and Fig. [S](b) shows the same for r = 
s — 2. The distributions are generated by eqn. Fig. [7] 
shows how the overlap distributions behave for s — r = 1 
when the ratio r/s goes to one. The overlap distribution 
with a randomized offset is generated by sampling over 
the overlap for 1000 different configurations and for all 
the possible offsets. 



V. SUMMARY AND DISCUSSION 

We model the static friction force between two atomically 
smooth surfaces by a two-chain version of Tomlinson's 
model. We consider substitutional defects in terms of 
vacancies along both chains (see Fig. [I}. The static 
friction force is assumed to be directly proportional to the 
number of beads in the upper chain which are locked in 
the potential wells in the potential arising from the chain 
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FIG. 5: (Color online) The collapse of the overlap 
distribution for different generations (n values), (a) for 
s, r = 3, 2 and (b) for s, r = 6, 5. The overlap 
distribution marked with lines are for random cantor 
sets, and are generated by eqn. The overlap 
distribution marked with points (marked by an asterisk 
* in the legend), are for random cantor sets with a 
randomized offset, and are generated by sampling. The 
x axis is in units of expected overlap. The embedded 
plot shows the same in log scale for the y-axis. 



beneath. This number is modeled by the overlap of two 
binary chains with a given correlation structure. We are 
in particular interested in the strong disorder limit where 
self-similarity may appear. The self-similarity translates 
to the power-law correlation in the spatial displacements 
of the remaining beads. The main motivation to consider 
this is the fact that the height profile of rough surfaces 
often has a fractal (self-similar) property (see e.g. p^j). 

A generalized random Cantor set is utilized to capture 
such properties. This Cantor set is generated by remov- 
ing all but r out of s segments from each remaining cle- 
ment at each generation. This procedure is explained in 
detail in section IIIII The remaining beads in the chain 
are distributed according to this random Cantor set. We 
study both the cases of with and without a random rela- 
tive offset between the Cantor sets (as illustrated in Fig. 
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(a) 




(b) 



FIG. 6: (Color online) The n invariant overlap 
distribution, generated by cqn. ((6]), for different values 
of s, r, with (a) s — r = 1, and (b) s — r = 2. the x axis 
is in units of expected overlap. 



The overlap of two chains is assumed to give the 
static friction force. Hence, the distribution of the static 
friction force is given by the distribution of the overlap 
of the Cantor sets. 

It may be noted here that the earlier applications of 
two fractal overlap models [l9H2l| . in the context of 
earthquake dynamics, focused on the time series of the 
overlap of regular (non-random) Cantor sets. Our in- 
terest here is the overlap of random Cantor sets in the 
context of static friction. 

Cantor sets have already been used to represent the 
scale invariance property of the contact area overlap be- 
tween two plastic surfaces of macroscopic objects in the 
study by Warren and Krajcinovic 22[. Their model docs 
not represent similar properties of each of the surfaces, 
as in our model. In the model presented here, Cantor 
sets (embedded in one dimension) are used to represent 
each of the surface and we calculate the overlap profile 
between them. Warren and Krajcinovic, on the other 
hand, use Cantor sets embedded in two dimensions, and 
calculate how it overlaps with a plane surface. As such, 



our model is different and the results are not comparable. 

For the overlap of the generalized random Cantor sets, 
without the randomized offset, we have found a recur- 
rence relation for the distribution (cqn. ((5])). Moreover 
we find that this distribution follows a scaling structure 
(cqn. Q). This scaling leads to a distribution which is 
independent of the generation n (though dependent on 
values of s and r). We further show how one can specify 
the distribution uniquely in terms of the moments, and 
do the calculation explicitly for the case when s = 3 and 
r = 2 (see eqn. ([5])). 

After introducing a randomized offset between the 
Cantor sets (with periodic boundary conditions), we no 
longer have a recurrence relation, but we find numeri- 
cally similar qualitative behavior. The distribution for 
the case with a randomized offset shows the same n in- 
dependent scaling behavior as the distribution without a 
randomized offset. Moreover, the distribution is shifted 
to slightly higher values of overlap as shown in Fig. [5] 
The same behavior is seen for different values of s and 
r. The embedded log plot in Fig. [5] shows that the tail 
behavior falls faster than exponential in both cases. 

When we look at the distribution for different values of 
s and r, generated by eqn. we see the emergence of 
multiple local maxima. This is a property of the sum of 
convolutions coming from cqn. ©, and depends on the 
allowed values of overlap for the first generation. These 
local maxima are averaged out when we look at the case 
with a randomized offset (see Fig. [S]and Fig. [7J). We have 
conveniently presented the overlap distributions in units 
of the expected overlap (the actual overlap of a given set 
is found by multiplying the x axis with \i n = (r 2 / s) n ). 

For a realistic chain for beads with vacancies, we can 
not assume that the remaining beads are distributed ac- 
cording to precise values of s and r. Nor can we assume 
that the overlap distribution is as for two Cantor sets 
without randomizing a relative offset. If we consider the 
limit where r/s approaches unity, (i.e. the limit where 
only a small fraction of beads are removed at every gen- 
eration), we find for the overlap distribution with a ran- 
domized offset, that the shape of the distributions have 
some common general properties. The distribution gets a 
peak at the expected overlap, but also an interval with a 
non-zero probability for values higher than the expected 
overlap, see Fig. 

The overlap of beads distributed as a Cantor set 
is qualitatively very different from the overlap of uni- 
formly displaced beads. To compare the distribution 
with that in Fig. [6j set p = (r/s) n and N = s" in 
cqn. ([TJ. The resulting distribution would then be ap- 
proximated by a Gaussian in the limit of large n, with 
mean Np 2 — (r 2 / s) n and a variance Np 2 (l — p 2 ) = 
(r 2 /s)"(l — (r 2 /s) n ). This would correspond to a sin- 
gle peak at x/fi n = I in Figs. [5][7J On the contrary the 
static friction force distribution for chain models, where 
the remaining beads are distributed in a scale invariant 
way, will not converge to a delta peak distribution but 
rather would be like that for Cantor sets with a random- 
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ized offset as in Figs. [5][7] 

It is hard to compare our results with experimental 
data on microscopic dry friction between surfaces hav- 
ing scale invariant disorder. Surfaces having microscopic 
self-affine disorder have been studied using AFM [23[ , but 
we are unfortunately not aware of any study where the 
static friction force between two such surfaces have been 
considered. We would like to mention that the above 
results for the distribution microscopic friction qualita- 



tively agrees with the observation of similar distributions 
of the friction coefficient of Aluminium alloys under cold 
rolling (see e.g., [24| ). Such an analysis can also be ef- 
fectively utilized for comparing the distributions for dry 
friction coefficients between plastic rock surfaces having 
well known scaling properties of asperities. However, in 
a different context of studying the effect of multiscale 
roughness on contact mechanics, similar analysis has al- 
ready been done (see e.g., p5|). 
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(a) 




FIG. 7: (Color online) The overlap distribution for 
s — r = 1 when r/s goes to 1. s, r = 7, 6 in (a), 8, 7 in 
(b) and 9,8 in (c). Lines (w/o randomized offset) are 
generated by eqn. ^j. Points (with randomized offset) 
are generated by sampling. The x axis is in units of 
expected overlap. 



